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DETECTION  AND  RETRIEVAL  OF  CIRRUS  CLOUD 
SYSTEMS  USING  AVHRR  DATA: 

VERIFICATION  BASED  ON  FIRE-II-IFO  COMPOSITE  MEASUREMENTS 

K.N.  Liou,  S.C.  Ou,  N.X.  Rao ,  and  Y.  Takano 
University  of  Utah 
Department  of  Meteorology/CARSS 
Salt  Lake  City,  Utah  84112 

ABSTRACT 

We  have  developed  a  detection  scheme  to  identify  single  and  multilayer  cirrus 
cloud  systems  based  on  the  physical  properties  of  the  AVHRR  Chs .  1-2  reflectances 
and  ratios,  the  brightness  temperature  differences  between  Chs .  4  and  5,  and  the 
4  brightness  temperatures.  Clear  pixels  are  first  separated  from  cloudy  pix( 
which  are  then  classified  into  three  types:  cirrus,  cirrus/low  cloud,  and 

cloud.  This  scheme  has  been  applied  to  the  NOAA  satellite  data  collected  over 
FIRE-II-IFO  area,  Kansas,  during  nine  overpasses  within  seven  observation  d 
(November  -  December  1991) .  We  have  validated  the  detection  results  against 
cloudy  conditions  inferred  from  the  collocated  and  coincident  ground-based  lidai 
radar  images,  balloon-borne  replicator  data,  and  NCAR- CLASS  humidity  sounding 
a  case-by-case  basis.  We  show  that  the  satellite  detection  results  are  consis 
with  the  cloudy  conditions  inferred  from  these  independent  and  complemer 
measurements.  We  have  also  modified  our  retrieval  scheme  for  the  determinati< 
cirrus  optical  depth  and  ice  crystal  size  in  multilayer  cirrus  cloud  systems 
case  study  using  FIRE-II-IFO  data  is  reported. 

1.  INTRODUCTION 

/ 

Cirrus  clouds  have  been  recognized  to  play  a  key  role  in  the  global  radiative 
energy  balance  and  climate  changes  (Liou  1986).  Information  on  cirrus  cloud 
parameters  is  critically  important  to  the  development  of  cirrus  cloud  formation 
models,  the  upgrade  of  real-time  global  cloud  analyses,  and  the  computation  of 
atmospheric  and  surface  radiative  parameters  in  climate  and  general  circulation 
models , 

In  recent  years,  our  research  group  has  developed  a  novel  and  comprehensive 
remote  sensing  algorithm  for  the  retrieval  of  cirrus  cloud  temperature,  optical 
depth,  and  mean  effective  ice  crystal  size  using  AVHRR  data  (Ou  et  al.  1993;  Rao  et 
al,  1995).  Validation  of  this  cirrus  remote  sensing  program  has  been  carried  out 
using  the  local  daytime  satellite  data  collected  during  FIRE-I-IFO  and  FIRE-II-IFO 
(Rao  et  al.  1995;  Ou  et  al.  1995a).  A  very  important  procedure  in  determining 
cirrus  cloud  parameters  is  the  detection  of  the  sky  condition  within  the  field- of - 
view  of  satellite  radiometers.  Our  detection  and  retrieval  schemes  have  been 
developed  primarily  for  applications  to  single-layer  cirrus  clouds. 

Surface  observations  show  that  multilayer  clouds  frequently  occur  in  the 
frontal  areas  where  cirrus  clouds  overlay  boundary  layer  convective  or  stratus 
clouds.  In  this  paper,  we  describe  a  numerical  scheme  for  detecting  multilayer 
cirrus  pixels  using  AVHRR  Chs.  1  (0.63  /xm) ,  2  (0.86  ^m)  ,  4  (10.9  |im) ,  and  5  (12.0 
/xm)  data.  Moreover,  we  also  present  a  preliminary  investigation  for  the  retrieval 
of  cirrus  cloud  optical  depth  and  ice  crystal  size  in  multilayer  cloudy  conditions. 
Verifications  of  the  detection  as  well  as  retrieval  schemes  utilize  the  composite 
data  sources  available  from  FIRE-II-IFO, 


2.  DETECTION  AND  RETRIEVAL  OF  MULTILAYER  CIRRUS  CONDITIONS 


2.1  DETECTION  differentiation 

^  -  with  the  availability  of  visible  ’can  be  made  using 

During  daytime,  -j^tions  over  various  types  of  must  be  less  than  a 

between  deer  end  cloudy  Ch.  1  ''Uectanc  ^ 

Che  following  ■=■^1'='’^^::.,!“  ''”'“el.eLnce“  f^ot  Che.  2  ani  1. 

threshold  value  *  ueces  y  ^cogtions .  to  identify 

we  can  define  a  tc'"-  ^  q  determined  from  the  ^  ^  Third,  the  Ch. 

.„d  use  «o  threaholds  ^.J.Xd'and  water  surf  aces  .respectivel^^^^^  ^ 

cloudy  and  clear  ^  clear  pixel  must  be  identification. 

4  brightness  a  -Id  temperature  can  be  establis  ^  functions  and 

&,-or‘ciUr^;  £»rand'. 

li:^tlly'Se":u:r“Uue.  they  are  1.-1  classified  .. 

after  all  the  “’■“"'‘y  cloud,  and  low  cloud,  pi,,cls  with 

three  classes:  cirrus  cltrus^/1^^  optically  thick  cirrus  the  visible- 

brightness  233  K  are  identified  as  ‘  than  those  for  cirrus 

temperatures  less  clouds  are  generally  Relatively  small  sizes 

channel  reflectances  for  cf  water  droplf"  ^he  latter 

clouds  because  the  “Jd  are  generally  ,  threLld  (-0.2)  can 

and  high  number  concen  reason,  a  visible-c  clouds.  Third,  the  Q- 

(Liou  1992,  Table  ^•2)-  pixels  that  contain  I®''  c^  clouds.  For 

be  established  to  filte  usually  smaller  than  the  .^.^us  over  low  clouds, 

ratio  for  low  clouds  larger  than  that  fo  cirrus  over  low 

Clttus  »-7^,tfo'Lr  citrus  over  Lter  Is  liw  cloud 

Moreover,  the  Q-ratlo  t  values  to  separate  cirr 

clouds.  Thus,  we  can  set  thresh  layer  cirrus  from 

or  clrrus/low  cl«ua.  c:,ed  to  e,eparate  dlngU^^j,^  .p,  .  threshold 

The  eioud  alone  conditions.  ^  5  pp  differentiate  the 

for  lo“  '^““^rse««te^rirWlow  cloud  and 


2.2  preliminary  retrieval  ,t-al  size  in  multilayer 

.  ,  .  the  cirrus  cloud  optical  depth  channels  follows  the 

,  <-v,p  aVHRR  0.63,  3.7,  and  .  /■1995)  and  Ou  et  al. 

cirrus  condition  using  the  1  (1993),  Rao  et  al.  retrieve  the 

:™erlcal  PP»P*r/%tr3  7  and  10  ’  -1  thermal  »dlances  are  used  to  retr^^ 

(1995a).  In  brief,  which  the  Ice  crystal  alaejnd^^^^ 

cloud  temperature  an  basis  of  cloud  mlcrophyslc  a  7  urn  radiance  for 

-".m\\er1nu.r  °a..oval  f 

r/oiSr/ 0.03  " 

iLgU  cirrus  dl«“1  „,re  couLted  during  nRE-II-I^^  of 

jrire  "rnfdmed  the  PP-^f  tSn  cirrfs^th^n  the  upwelUng  radlancej 

low  cirud.  If  Ita 

reaching  the  sinele-layer  cirrus  case.  r-f^moval- retrieval  progr 

rtriner/tr  Se^^v  sS  tadl^^j;;  Thotr  cS^ur^iow  clonds  have  th, 
developed  by  Rao  et  al.  (19y:5:>- 


P  coverage,  information  of  the  thermal  upwelling  radiances  in  the  3.7  and  10  9 
^  hannels  as  well  as  the  low  cloud  albedo  is  unknown  and  must  be  assumed  a  priori. 

rhis  case  we  use  the  climatological  microphysics  data  for  stratus  to  perform 
Iheoritical  calculations  to  obtain  the  required  inputs  in  retrieving  the  cirrus 
I^tical  depth  and  ice  crystal  size.  We  can  then  compute  the  visible  radiances  at 
top  of  the  atmosphere  and  compare  with  observed  radiances  to  assess  the 
^liabUity  of  the  calculated  optical  depth  for  low  cloud.  Subsequently  iterations 
be  devLoped  to  derive  a  consistent  set  of  optical  depths  for  both  cirrus  and 

]^ow  clouds . 

3  VALIDATION  OF  THE  DETECTION  AND  RETRIEVAL  SCHEMES  USING  FIRE-II-IFO  DATA 

The  FIRE-II-IFO  was  carried  out  near  Coffeyville,  Kansas,  during  November  and 
n.rember  1991  There  were  a  number  of  dates  during  which  multilayer  cloudy 

For  the  detection  o£  cirrus  cloud  plKels,  the  high- resolution 
fHRPT)  AVHRR  data  from  NOAA-11  and  NOAA-12  polar-orbiting  satellites  are  used^  We 
Lve  acquired  ground-based  lidar  and  radar  images,  balloon-borne  replicator  data 
and  NCAR-CLASS  humidity  soundings  on  a  case-by-case  basis.  From  ne 
datasets  we  have  selected  seven  representative  dates  (nine  overpasses)  for 
study  including  clear,  cirrus,  and  cirrus/low  cloud  conditions.  For  each  case  we 
compare  cloud  tyves  identified  from  satellite  radiances  with  those  derived  from 

ground-based  composite  instruments.  o  t  noi-  i-Vie 

^  Table  1  summarizes  the  results  of  the  comparisons.  Columns  2-6  list  the 

required  parameters  for  numerical  processing  determined  from  satellite  data,  while 
columns  7^-11  depict  cloudy  conditions  obtained  from  satellite  data  and  various 
gri““bas.d  and  in  situ  Inatrcanta.  Overall,  the 

fuccessfully  differentiates  among  clear  (12/6b) ,  cirrus  (ll/26b  and  12/5  ) , 
cirrus/low  cloud  conditions.  These  results  are  consistent  with  the  clou  y^ 
conditions  identified  from  the  independent  and  complementary  ground-based. 

measurements^  selected  the  cirrus/low  case  that  occurred  on  29  November  1991  for 
testing  the  retrieval  scheme.  Figure  1(a) 

humidity  profiles  obtained  from  the  NCAR-CLASS  sounding  launched  at  1343  UTC.  The 
Sius  cloud  base  and  top  heights  derived  from  the  replicator,  PSU  94  GHz  cloud 
radar,  and  visible  lidar  are  -  6(-20‘>  C)  and  9  km  (-41°  C)  respectively  A  moist 
layer  roughly  corresponding  to  the  cirrus  cloud  layer  is  evident.  Moreover  an 
moist  layer  (RH  >  90  X)  existed  between  -  1  and  2  km,  corresponding  to  ^ 
cloud  layer  detected  by  the  PSU  radar.  Temperature  inversion  occurred  at  ^e  low- 
level  clLd  top  and  near  the  peak  of  the  relative  humidity  around  7  km.  ^ 

retrieved  cloud  temperature  is  233  K,  which  is  the  average  of  62  pixels  within  the 
0  1°  X  0  1°  area.  The  standard  deviation  is  5.8  K,  indicating  that  clou 
temperatures  were  not  uniform  within  the  retrieval  domain.  Moreover,  the  mean  cloud 
height  determined  from  the  temperature  sounding  is  9.2  km,  which  is  near  the  cloud 

On  29  November  1991,  there  were  only  three  levels  of  replicator  measurements 
available.  At  9.13  km,  ice  particles  are  composed  of 

irregular  crystals  with  the  maximum  dimensions  ranging  from  25  to  425  (xm.  The  size 
distribution  peaks  at  75  /xm  with  a  number  concentration  of  0.24  L  M"!  ^  At  8_18 
km  the  ice  crystal  size  distribution  is  similar  to  that  at  9.13  km.  However,  the 
l^par  lliit  ./th,  measured  sizes  Increases  to  575  pm.  The  size  distribution  also 
peaks  at  75  /xm.  but  with  a  smaller  value  of  number  concentration  of  0.1  L  /xm  . 
The  level  at  7.46  km  corresponds  to  a  local  peak  in  the  relative  humidity  pro  i  e. 
The  size  distribution  with  a  upper  limit  of  875  /xm  is  broader  than  the  previous  two. 
The  ice  crystal  shapes  include  bullet  rosettes  and  aggregates  The  derived  mean 
effective  sizes  for  these  three  levels  are  96,  116,  and  146  /xm  from  top  to  bottom 
The  vertically  averaged  mean  effective  ice  crystal  size  is  134.6  /xm  (solid  vertical 
bar  in  Fig.  lb)  .  The  retrieved  mean  effective  size  from  satellite  radiances  is 
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1,  M,.  rpolicator  derived  and  the  satellite 

130  9  /im.  On  the  bottom  scale  are  ®  respectively. 

retrieved  optical  depths,  which  are  2_-21  tL  retrieval  program  for  the 

„e  are  Spth'nd  n,.a„  ice  cryat.l  aUe.  ««re  eo»prehensrve 

iSy"eranrv.ll<iatlot,a  will  be  reported  In  the  future. 
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Table  1 . 


Results  of  the  satellite-based  cloud  detection  compared  with  ground-based  radar,  lidar, 
and  balloon-borne  replicator  measurements. 


r,:  Ch.  1  reflectance 

Q:  Ratio  of  Ch.  2  to  Ch.  1  radiances 

BTD45:  Brightness  temperature  difference  between  Ch.  4  and  Ch.  5 
T4:  Ch.  4  brightness  temperature 
a:  NO  A  A- 12  overpass  (-  1400  UTC) 
b:  NOAA-II  overpass  2100  UTC) 


Parameter  Values  Cloudy  Condition 


Date 

number 
of  pixels 

r.(%) 

Q 

BTD45(K) 

T,(K) 

satellite 

PSU 

radar 

ETL 

lidar* 

LaRc  soundings 
lidar* 

12/6b 

117 

12.1 

1.22 

0.92 

287.0 

clear 

clear 

clear 

clear 

— /dry 

12/5b 

89 

32.1 

1.07 

3.04 

249.4 

cirrus 

cirrus 

cirrus 

cirrus 

ice/dry 

11 /26b 

79 

24.2 

1.10 

2.73 

271.6 

cirrus 

cirrus 

cirrus 

cirrus 

ice/dry 

11 /22a 

93 

57.7 

0.91 

0.46 

244.7 

ci/low 

ci/low 

cirrus 

cirrus 

ice/low 

11 /29a 

52 

45.6 

0.89 

1.18 

249.7 

ci/low 

ci/low 

— 

— 

ice/low 

11 /28a 

100 

23.8 

0.91 

2.09 

272.1 

ci/low 

ci/low 

cirrus 

cirrus 

ice^/low 

11 /28b 

85 

20.0 

1.04 

1.77 

284.2 

cirrus 

•  cirrus 

cirrus 

cirrus 

iceVdry 

11 /27a® 

‘  5200 

44.5 

0.91 

0.80 

246.0 

ci/low 

ci/low 

— /low 

11 /27b 

67 

63.5 

0.93 

3.74 

262.7 

ci/low 

ci/low 

— /low 

*•  Both  ETL  lidar  and  LaRc  lidar  measured  signals  from  the  backscattering  of  boundary  layer  aerosols  and 
low  cloud  particles.  These  signals  were  not  included  in  ^e  images  analyzed  in  this  study. 

*f '  Based  on  replicator  measurements  between  the  two  satellite  overpasses. 

Satellite  cloud  detection  results  are  based  on  data  over  l.CPxl.CT  area  around  SS.S**  N,  96.5^  W. 
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Figure  1.  (a)  Cloud  ba«c  and  top  heighu  for  cirrus  and  stratus  determined  from  lidar,  radar  and  sounding  daU,  as  well  as  temperature  and 
humidity  profiles  obuined  from  the  NCAR  CLASS  sounding  system  at  1343  UTC,  29  November  1991,  Overlapped  with  the 
temperature  profile  are  the  mean  retrieved  cirrus  cloud  temperature  over  O.l*  x  0,1®  domain  around  Coffeyville,  Kansas,  and  (b) 
Di^lay  of  the  itipUcator-dcrived  mean  effective  sizes  at  selected  height  levels,  their  vertical  average,  and  the  retrieved  value.  Also 
shown  on  the  bottom  scale  are  the  <^tical  depths  derived  form  the  replicator  data  and  from  the  retrieval. 
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